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We consider the van der Waals interaction between a pair of polymers with quenched hetero¬ 
geneous sequences of local polarizabilities along their backbones, and study the effective pairwise 
interaction energy for both stiff polymers and flexible Gaussian coils. In particular, we focus on the 
cases where the pair of polarizability sequences are (i) distinct and (ii) identical. We find that the 
pairwise interaction energies of distinct and identical Gaussian coils are both isotropic and exhibit 
the same decay behavior for separations larger than their gyration radius, in contradistinction to 
the orientationally anisotropic interaction energies of distinct and identical stiff polymers. For both 
Gaussian coils and stiff polymers, the attractive interaction between identical polymers is enhanced 
if the polarizability sequence is more heterogeneous. 

PACS numbers: 05.20.-y, 66.30.hk 


I. INTRODUCTION 

Van der Waals (vdw) forces [THll are long ranged forces which act between all kinds of bodies including electrically 
neutral ones, and are therefore prevalent in nature. These forces lead to a variety of phenomena including flocculation 
in colloid systems, [S] and are also responsible for the ability of geckos to stick to walls. [S] Van der Waals forces 
furthermore appear to be active in the first stages of planet formation, in causing planetisimals to cohere together 
when the gravitational force is still too weak to bring about cohesion. [7] In industry and water treatment plants, vdw 
forces have also been utilized to treat waste water. Waste water contains impurities which are charged, and coagulants 
are added to make them neutral. These neutral impurities which are colloidal then attract each other by vdw forces 
and settle down as Hoc which can then be easily disposed of. [8] 

More relevantly for evolutionary biology and the physics of self-assembly, vdw forces also contribute to the process 
of recognition between molecules. Molecular recognition is the process whereby a (bio-)molecule, say a protein, 
recognizes a specific structural feature on another macromolecule that it interacts strongly with. [giin! This includes 
the recognition between proteins and DNA molecules that are so essential for biological functions. A striking feature 
of biological molecular recognition is its high target sequence specificity, [nmu for example, certain protein molecules 
such as the lac repressor are able to bind to a specific sequence out of six million possible ones on the DNA of an E. coli 
bacterium. On the other hand, in molecular electronics and DNA computation, [13] there has been increasing interest 
in the possibility of self-assembling DNA-based circuit boards, the self-assembly being facilitated by the molecular 
recognition between pairs of DNA molecules. m An investigation of vdW interactions between polymers would thus 
be of relevance to better understanding the process of molecular recognition. 

Our paper is structured as follows. We first look (heuristically) at how vdW interactions can arise between neutral 
but polarizable atoms. We then consider physical aspects of the dsDNA molecule, which motivates our study of 
vdW interactions between stiff polarizable polymers. Next, we explain the formalism we use to describe such vdW 
interactions, focussing on the features that can aid the mechanism of molecular recognition. m Towards the end we 
will also consider the case of flexible polarizable Gaussian coils (which can describe globular proteins such as lysozyme 
and RNA polymerase), and how their vdW interaction differs from that between stiff polymers. 


A. Van der Waals interactions: Heuristics 

First, let us consider heuristically how vdW interactions can arise. Let us consider a neutral atom, whose proton 
charge balances the charge of its electron cloud. At finite temperature, the electron cloud undergoes thermally driven 
distortions, and at any one instant the center of the electron cloud will be displaced relative to the position of the 
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proton. It thus acquires an induced dipole moment. A second atom in the vicinity of the first atom experiences an 
electric field generated by the induced dipole of the first atom, which polarizes the electron cloud of the second atom, 
giving rise to another induced dipole. Although each induced dipole vanishes under time averaging, the product of 
two induced dipoles does not vanish, and this gives rise to a non-zero polarization-type or van der Waals interaction 
between the atoms. 


B. Structure of dsDNA molecules 

Next, let us look at some structural features of dsDNA molecules that are relevant to molecular recognition. [16] 
The molecules have persistence lengths of around 100 base pairs (about 50 nm), which means that short dsDNA 
molecules (e.g., 10 to 20 base pairs long) essentially behave as stiff rodlike molecules. Each dsDNA molecule, having 
a length that is much greater than its radius, is correspondingly much more polarizable along its backbone than 
in the directions transverse to it. Additionally, each dsDNA molecule comes equipped with a sequence of different 
nucleotide base pairs, different base pairs having different polarizabilities depending on the type of base pair and 
the identity of its neighbours. In solution the highly charged dsDNA molecule also attracts counterions to the 
backbone, which further modifies the local polarizability along the backbone. All this implies that the associated 
sequence of polarizabilities is heterogeneous along the backbone of the dsDNA molecule. Significantly, as the base 
pair sequence is structurally built into the backbone of the molecule, the corresponding sequence of polarizabilities is 
also thereby quenched (i.e., the polarizabilities are effectively hxed, albeit random, on experimental time-scales). In 
general there will also be an annealed (i.e., thermally fluctuating) contribution to the local polarizability arising from 
the condensed counterions, but for simplicity, we only consider the situation in which the polarizabilities are entirely 
quenched. 

The base pairs are hydrophobic, which means that they carry induced dipoles and do not carry permanent dipoles. 
On the other hand, the sugar phosphate backbone of the dsDNA molecule is highly negatively charged with a linear 
charge density of around 6 elementary charges per nm. In solution, this induces a cloud of counterions to condense 
on the backbone bringing down the net charge of the dsDNA to an effective charge of around 1.4 elementary charges 
per nm in the case of monovalent counterions. [Mns] The electrostatic potential due to this effective charge is further 
screened through the Debye screening mechanism. Associated with the latter effect is the Debye screening length, 
defined by Ad = {££okBT/J2j ^^9^), where e is the relative permittivity of the solvent medium (for water, e « 80), Eq 
is the vacuum permittivity, is Boltzmann’s constant, T is the temperature, Cj is the concentration of ion species j, 
and Qj the corresponding charge valence of the ion species. For example, at room temperature (T = 298 K) and for a 
salt solution containing only NaCl of concentration IM the corresponding Debye length is given by Ad ~ 0.3 nm. At 
separation distances greater than Ad, electrostatic interactions are much less important than vdW interactions, and 
the molecule is effectively neutral. We consider inter-molecule separations of order 1 to tens of nm, which means that 
we can essentially disregard electrostatic interactions and focus only on vdW interactions. 


C. Objective and strategy 

In what follows, we are going to address the following questions: (i) How does the vdW interaction between two 
polymers depend on the correlatedness of their polarizability sequences? Specifically, is the vdW interaction statisti¬ 
cally stronger between a pair of polymers with identical polarizability sequences than between a pair of polymers with 
distinct sequences? (ii) What effect does a more pronounced polarizability heterogeneity have on vdW attraction? 
(iii) How does the vdW interaction depend on the relative orientation of two polymers? We will see that the vdW 
interaction between a pair of stiff polarizable polymers is characterized by the interplay between polarizability se¬ 
quence specificity, heterogeneity, and orientational anisotropy, all of which contribute to the mechanism of molecular 
recognition. 

For the rest of the paper, we adopt the following strategy. Firstly, we approximate a dsDNA molecule by the 
shish-kebab model. |2D] In this model, the polymer is a stiff rod, and each monomer (i.e., a nucleotide base pair) 
is assumed to have the same size a. The monomers are represented in Fig. [^ by circles with different shades, 
corresponding to differences in polarizability. Each monomer can have a different polarizability value. As dsDNA 
molecules are non-zwitterionic, we only consider interactions between induced dipoles. We have also approximated 
the polarizability by its static value, an approximation that works well at high temperature. On the other hand, we 
consider only non-retarded interactions, i.e., corrections due to the finite speed of light can be neglected, and this is 
a good approximation at inter-molecule separations smaller than the retardation lengthscale fic/feeT (obtained, e.g., 
by comparing the London dispersion formula and Casimir-Polder formula for a pair of atoms). Thus we see that the 
non-retarded approximation is only valid for temperatures that are not too high. For biological systems, T « 300 K, 
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FIG. 1: “Shish-kebab” model of a dsDNA molecule as a stiff polymer consisting of monomers (i.e., nucleotide base pairs, 
represented by circles) of equal size a. The differences in shade reflect differences in polarizability. 


which corresponds to a retardation lengthscale of 7 /im. The range of separations that we are interested in probing 
(one to tens of nm) thus falls safely within the non-retarded regime. It is well known that water has structure (water 
has also been regarded as a network of hydrogen bonds). To make the problem analytically tractable, we approximate 
the aqueous medium as a homogeneous dielectric of e « 80. 


II. HAMILTONIAN FOR INDUCED DIPOLES 


A. Dipole electrostatics 


To describe the vdW interaction between polarizable polymers, we proceed by constructing a Hamiltonian for the 
induced dipoles. This requires us to first define the dipole charge density. In a neutral molecule, if the charge of the 
proton is q and the separation vector between the proton and the center of the electron cloud is x, the dipole moment 
is given by p = qx. The corresponding dipole charge density is given by p(r) = —qS{r — R + |) + qS{r — R — |), 
where R denotes the position vector of the centre of the dipole. For small d, we can do a Taylor expansion. This 
yields p{r) « —p • V(5(r — R). For N polymers each with M monomers, the dipole charge density is given by 

N M 

pW = “X! X! ■ vd(r- RW(m)), (2.1) 

i—\ m—1 

where now p^*^ (R*^*^(to)) refers to the induced dipole moment (position vector) of monomer m on polymer i. Having 
defined the dipole charge density, we can write down the standard expression for the potential at position r, viz., 


$(r) 


1 

dweso 


;3r P(r' 

|r - r'l 


( 2 . 2 ) 


The electrostatic energy due to dipoles is given by 

ffp,N = d^rp{r)^r) 

3(p(^)(m) ■Rij-(TO,n))(p(j)(n) ■R,j-(m,n)) ) 

(i?,,(m,n))5 /• ^ 

Here Rij(m,n) = R^-^)(n) - R*'®^(to), the factor of 1/2 is to prevent double counting in the sum over pairwise 
interactions of dipoles. By plugging in our formula for the dipole charge density, we find that the interaction potential 
between two dipoles is equal to two gradients acting on the Coulomb potential between two monopolar charges. 


B. Dipole “elastics” 

In addition to the foregoing electrostatic term describing the interactions between dipoles, we also have an elastic 
term describing the energetic cost of distorting the electron cloud of each molecule. If we have an external field E, 
the field will induce a dipole moment p, and these quantities are related via the equation p = a • E, where a is the 
polarizability tensor. A larger value of a means that the molecule is more polarizable. Rewriting E = a~^ ■ p, we 
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see that this relation is reminiscent of Hooke’s law, and thus we can regard a~^ as a “spring constant” and p as 
an “extension.” As the polymer is much more polarizable along its backbone than in directions perpendicular to it, 
we can approximate E Ri where a is now the polarizability along the backbone. By analogy with Hookean 

elasticity, we can write the “dipole elastic distortion energy” for a monomer, viz., 

ffe = (2.4) 

For N polymers each with M monomers, the dipole elastic distortion energy is given by 


N M 


He,N 


ah) (to) 


2=1 TYl — l 

The full Hamiltonian for N polymers is then given by 

Hn = 


(2.5) 


( 2 . 6 ) 


As biological systems are typically at room temperature, there are many different thermal configurations of the induced 
dipoles, each characterized by a different Boltzmann weight. The effective interaction energy between a pair of rigid 
(i.e., conformationally non-fluctuating) polymers comes about by coarse-graining over their respective induced dipole 
degrees of freedom. As the dipoles are much more easily induced along the polymer backbone, we can equate the 
direction of the dipole to the orientation of the polymer, given by the tangent vector For stiff polymers, the 
direction of each dipole is constrained by the orientation of the polymer, and thus for a given polymer conformation 
only the magnitudes of the dipoles undergo thermal fluctuations. Thus we can write the partition function for a pair 
of rigid stiff polymers as 


Z(iV,{R(*)}) = n / dp^"Hrn)e-<^^^. (2.7) 

2,m 

It is very difficult to evaluate the partition function exactly, owing to the large number of interacting dipoles. One 
method of making the problem analytically tractable is the Hubbard-Stratonovich transformation method (see, e.g.. 
Ref. [21]), where we introduce a local auxiliary field ip to decouple the dipoles. The dipoles can then be integrated 
out, and what remains is a local theory for (j9.|15j viz., 


Z(fV,{RW}) cx 




— j/3 X] ("*))) 


( 2 . 8 ) 


The resulting Hamiltonian is quadratic in the gradient of p. We can interpret p as the local field fluctuation generated 
by the presence of the induced dipoles. We also see that the dipoles effectively contribute to the dielectric permittivity 
of the polymer: away from the polymer, the dielectric permittivity is eeo, and at the position of the polymer, there 
is a correction proportional to the dipole polarizability. 


III. EFFECTIVE INTERACTION ENERGY 


The calculation of the partition function can be simplified if the polarizability is small. We can then expand the 
exponent in powers of a up to second order, the second order being required to see the interaction between the pair 
of polymers. This is also formally equivalent to making a Hamaker summation approximation. We find|15) 

Z(iV,{RW}) « (3.1) 


where is the sum of the free energies of individual polymers and Tint is the effective interaction energy between 
pairs of polymers. For a pair of polymers, we find 

f F(l)(si) • V^^Hs2) 3(u(l)(si) • Ri2)(u(^Hs2) • R 12 ) F TO 

I Rh J ’ ^ 
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where R 12 = r — and are the tangent vectors of the first and second polymers, r is the 

separation between their centers of mass, and si and S 2 are the arc lengths of the first and second polymers. In 
deriving the above result, we have taken the continuum limit, where the sum over segments becomes an integral 
over the arc length. We note that inside the integrand the interaction decays as which is the usual decay law 
associated with non-retarded dispersion interactions. 

In the effective interaction energy, we also need to average over the polarizabilities. We assume that the 
local polarizabilities follow a Gaussian distribution, and the local polarizability is composed of a mean and a 
quenched fluctuation, i.e., = ao + The mean ao and variance specified by ((a^*^(s))) = ao and 

(s)(5a^'l^ (s'))) = g^aS{s — s')Sij, where ((...)) denotes a sequence average over the statistics of the polarizability. 
From such a statistics, we can deduce the following two limiting cases. For a pair of polymers with distinct sequences 
(i.e., a(*^(s) ^ a^^\s)), this implies that ((a^* )(s)aW(s'))) = Og, whereas for a pair of polymers with identical se¬ 
quences (i.e., a W(,) = a^^\s) = q;(s)), the polarizability correlator receives an additional contribution from the 
variance, viz., lla^'‘\s)a^^\s')'f} = ag + g'^a5{s — s'). 

In what follows, we first give an overview of the main results of Ref. m, for the interaction behavior of two stiff 
polymers, in the two limiting cases where the polymers have distinct and identical polarizability sequences. We then 
consider the interaction behavior of two flexible polarizable Gaussian coils. 


IV. RESULTS 


A. Two stiff polymers 


Let us first consider stiff polymers with distinct sequences. In the near-held regime (by which we mean their 
separation distance is much smaller than the length of each polymer), we can approximate each polymer by an 
inhnitely long polymer. From Eq. (3.2), we have found that the sequence averaged interaction energy is given hym 


{{Fint)) « - 




64^(eeo)2^2 vT^ m 


(4.1) 


where z = £ is the length of each polymer, and R );2 is the shortest length vector that connects the two 

polymers. In the far-held regime (i.e., where the separation distance is much greater than the length of each polymer), 
we have found that [15] 


mnt)) « - 


(z - 3yij/2)^M^Q^fcB'T 
327r2(£eg)2r6 


(4.2) 


where yi = ■ xjr and ?/2 = ■ r/r. Note that in both the near- and far-held regimes, ((Tint)) is invariant with 

respect to inversion of each polymer about its center, as there is no correlation of polarizability huctuations for the 
two distinct polymers, and thus each polymer sees a uniform (sequence averaged) polarizability on the other polymer. 

Next, we consider the sequence averaged effective interaction energy for a pair of identical sequences. As before, 
we hrst look at the behavior in the near-held regime. The sequence averaged energy is now given by 


((Fint)) -Fo + SF, 


where 


Fn = - 


6F = - 


M^alksT 

647r(£eo)2£2 7fr^|R*2|4> 

3kBTMg‘^xiyi,y2,z) 

16384•\/27r(££g)2fr® 


(4.3) 

(4.4) 

(4.5) 


We hnd two contributions: the hrst (Fg) is the contribution coming from the mean polarizability, and is identical with 
the sequence averaged interaction energy of two polymers with distinct sequences, whereas the second contribution 
(SF) is new, and comes from the correlation of quenched polarizability huctuations. Note that 5F decays as r~^. In 
the result above, y describes the orientational anistropy of the vdW interaction, and is given by 


xiyi,y2,z) = 


yr^(i-7)"^' 

5(3 -I- lOz + 3z^)(^i - 1 -^ 2 )^ 
1-7 


9 -b 14z + 41 z2 

105(yi +y2)^ 


+ ■ 


4(1-7)" 


(4.6) 
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Here, 7 = distinguishes between polymers that are aligned and those that are anti-aligned. We can 

distinguish between a state of alignment and a state of anti-alignment, because now the sequences are perfectly 
correlated, and to flip one sequence you have to flip the other too to get the same energy. 

In the far-held regime, we hnd that the sequence averaged effective interaction energy is given by 


{{Pint)) 


+ Mg'^){z -iyiy2Y 
327r^(e£o)^r'® 




(4.7) 


where we see that the only effect of the polarizability correlation is to renormalize the overall prefactor. 


B. Two flexible Gaussian coils 

We now address the question of how the interaction behavior is modihed when the pair of rigid stiff polymers are 
replaced by a pair of hexible Gaussian chains. For such chains, the local tangent vector of the coil and the 

polymer coordinate are independent of each other, in contradistinction to the case of stiff polymers. Let us 

write Fint for the effective interaction energy for a given chain conformation. Recognizing the complexity of a full- 
hedged description of interaction between given conformations of Gaussian chains, the problem can be simplihed by 
going to a coarse-grained description of interaction between Gaussian coils, where the thermal averaging over polymer 
coordinates and local tangent vectors has already been carried out on the partition function Z. This leads to 


= Yl /dRW(s)/dFW(s)Z[{R«(s),F«(s)}] 


OC 1 /3 ( (itself ){R(0} + (^int){R(0}) 


(4.8) 


In the second step we made use of the second order approximation as in Eq. dni ), and in the third step we re¬ 
exponentiated the terms. The overhead bar denotes averaging with respect to the local tangent vectors, and (.. -({rco} 
denotes averaging with respect to polymer coordinates The sequence average can also be performed on 

(Fint){R(0}- 

In the high temperature regime, we can make a further simplihcation by making use of the rotating dipole 
approximation, [55] which consists of (isotropically) averaging Tint over all orientations in space. This approxima¬ 
tion is valid when the Gaussian coils are spherical, which is the case when the centers of mass of the coils are 
sufficiently far apart, i.e., by a distance larger than the effective radius of gyration of the coils (defined in Eq. (4.20)). 
To perform the orientation average, we make use of the following relations valid for Gaussian chains: 


ta\s)t^^\s') = 0 ; 


where a,b = 1,2,3 are Gartesian indices. After averaging Fi„t [cf. Eq. over the orientations, we obtain 

3fcBT f ds f ds'a^^\s) 


Fint — — 


167r2(e£o)2 J a 


iRl2{s,s'W ’ 


(4.9) 

(4.10) 

(4.11) 


where Ri 2 (s,s') = R^^^(s') — R^^^(s). We still need to perform the average over fluctuating polymer coordinates 
as well as the sequence average over fluctuating polarizabilities. Parallel to our analysis for the stiff polymers, we 
can study the interaction behavior of a pair of Gaussian coils when (i) their separation distance is much larger than 
Rg (the effective radius of gyration of the coils, which is of nanometer order, cf. Eq. (4.20)) but still smaller than 
the retardation lengthscale (which is of micron order at T = 300 K), and (ii) when the coils are near each other at 
separations > Rg, in order that the Gaussian coil conformations still follow an isotropic Gaussian distribution. For 
each regime, we consider the following two cases: (i) distinct and (ii) identical polymers. 


1. Coils far apart 

We consider the regime in which two coils are separated by a distance much greater than their effective radius of 
gyration. Let us first study the case of distinct polymers, for which there is no correlation of polarizabilities between 
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the coils. The sequence average of the interaction free energy in Eq. (4.11) thus yields 

' ds f ds' 


{(F-i)) = - 


1 


167r2(£eo)^ J o-J o (^ 12 (S;S'))® 


(4.12) 


The centers of mass of the first and second Gaussian coils are given by and Rg^ = f ^R*^^^(s). Let 

us define the fluctuations JR(*)(s) = R*^*^(s) — R^^ and (5R(s, s') = Ri 2 (s, s') — (Rq^ — Rq^) = i5R^^^(s') — 5R^^)(s). 
In the far-held regime, SR <SC d, where d = r[?^ — R^^ is the separation between the centers of mass of the Gaussian 
coils, and we can perform a Taylor expansion (to quadratic order in ^r): 


{{Fint}) 


SkBTal fds fds' f 3|.5Rp 24(d ■ JR)^ ) 

167 r2(£eo)^d® J a J a \ # J 


(4.13) 


To make further progress, we assume the thermal huctuations dR^^^ and hR^^^ are Gaussian distributed, and average 
{{Fint)) over these huctuations. The radii of gyration of the individual coils are given by and which are 

dehned by (Rg^)^ = /^(|5R^*^(s)p)ii;.. For isotropic coils, we have 


I j{SRi^\s)SRi^\s))n, = i< 5 .,(RW)^; 

J ^(<5rW(s')<5rP^(s'))«, = i<5a.(Rf )^; 

I j{dR<~^\s))n, = {SR^^Hs'))n,=0. 


Performing the average over polymer coordinates, the interaction free energy becomes 


(((Rint)))fli,fl2 = - 


BfcB'TggM^ 

IGtt^ (££ 0 ) 2^6 


1-f 5 


(i?W )2 + (i?( 2 )) 2 ' 


d 2 


(4.14a) 

(4.14b) 

(4.14c) 


(4.15) 


The vdW interaction between a pair of Gaussian coils is thus attractive (like in the case of stiff polymers), but in 
contradistinction to the stiff polymer case the interaction is isotropic, which we expect because of the sphericity of 
the coils. 

Now let us consider identical Gaussian coils. The polarizabilities are now correlated, and the variance of polariz¬ 
ability huctuations generates an extra contribution to the sequence average of Eq. (4.11[). Gonsequently, we hnd 


(((^int))) 


Rl ,i?2 


3kBT{M^al + Mg^ 
167r2 (££ 0 ) 2^6 


1 + 5 


(4^^)^ + (Rf)2 

d 2 


(4.16) 


Similar to the case of identical stih polymers, the effect of the identicality of the Gaussian coils is an overall renor¬ 
malization of the ehective polarizability. 


2. Coils near each other 


It is non-trivial to study the interaction of two Gaussian coils whose separation is smaller than the effective radius 
of gyration, owing to the asphericity of the coils. On the other hand, it is possible to estimate the effective interaction 
energy of two coils that have a separation distance comparable to their effective radius of gyration, under the condition 


that d > 4i?g. To proceed, let us rewrite the dipole-dipole interaction kernel from Eq. (4.12) as follows: 


’ ds 
a . 
' ds 
a 

' ds 
a . 


'ds' 


a |R(2 )(sO-R(iHs)I® 

i3 


'+/d-r 

a 


3 <I(r-RW(s))<5(r-R(2)(+)) 

|r' — r |6 


[d^r 

a J J 




d-^r' 


d?k f dr’k' 


{2tt)^ J {2tt)^ 


(4.17) 





















In the second step, we have used the integral representation of the Dirac delta-function. For small fluctuations (5R^*^(s) 
we can expand the exponent to quadratic order in the fluctuations. Carrying out the thermal average over (s) 
and re-exponentiating, Eq. (4.17) becomes 


ds f ds' 


— — d^r d-^r 




d^k f d^k' ) 


(27r)3 J {2tt)^ 


1 - 


' ds 


{SRi^Hs)SRl^\s))n, - / ^{SR^^\s')SR<f\s')) 


R2 
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d^r 


d^r' I 


(27r)3y (27r)3 |r'— r|® 


Slr-R^b^ 3|r'-R,^h^ 


27M2 


8^3(4i))3(42))3 


d^r / d?r' ■ 




2 (Hrb 


(2)s2 


r' - r' 


(4.18) 


The interaction is between two polymer segments at r and r', belonging respectively to the first and second Gaussian 
coils, and the interactions are weighted by the Gaussian probabilities of finding the corresponding polymer segments 
away from the centers of mass of their coils. The integration is facilitated by making the following change of variables: 
r = c — f A, r' = c -I- f A. We have 


’ ds f ds' 


1 


a \|RG)(sO-RW(s)|6/i^„fi, 


27M^ 


d'^c 


'd^X - 


3|c-iA-R,' '|2 3|c+ii-R^'|2 






87r3(i?«)3(i?f )3 J 

27M^ /27ry/^ (R^^^)3(R^^))3 7 ^3A - 

. 3 J 3((4^^)2 + J 


87r3(i?«)3(42) )3 


A-d|- 

2-R2 


m 2 


M3A 


2h2 


247r3/2i?3/ A6 


(4.19) 


In the above, we have defined the effective radius of gyration Rg of the two coils, viz., 

Rl=\{{R^PY + {Rf^f). 


(4.20) 


At room temperature, the radius of gyration is typically of nanometer order for globular proteins; for example, the 
radius of gyration of lysozyme is approximately 2 nm and that of RNA polymerase is approximately 5 nm. [23] In the 
second step, we have integrated over c. We see that the integral over A is essentially the dipole-dipole interaction 
averaged over a Gaussian probability d istrib ution of sep aration distances that is peaked at a value A = |R[;^ — R[j^ |, 
with a variance 5\^ = R?. Using Eqs. (4.12) and (4.19), the interaction energy can be simplified to 


— k^TalM^ 

AFint)) — 


-d^X -A- 


I287r^/2(£eQ)2i73 J \6 


(4.21) 


Let us rescale A —>■ i?„ A, so that A is dimensionless. We can write 


'd^X -J 


R 


-3 


'd^X iA-d |2 

2 =2Tr Rg 


dx 


dXe y 


(4.22) 


where we have also rescaled d —)■ di?g (so that d is now dimensionless); x is the cosine of the angle between A and d, 
and / is given by 


f = ilnX + ^{X^ + d^ - 2xXd). 


(4.23) 


Next, we make a saddle point approximation to JdxJdXe We approximate the integral by replacing it with 

the configurations that have / minimized with respect to x, i.e., x = 1. Of these configurations, we approximate the 
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FIG. 2: Behavior of e~^ (cf. Eq. (4.231) as a function of \/Rg for the values |Rq^ — = 5,x = 1. We note the existence 

of two stationary points: a local maximum (corresponding to A+) and a local minimum (corresponding to A-). If the centers 
of mass of the two Gaussian coils are farther than Rg, the dominant contribution to will come from A+. 


integral over A by the solution to the saddle-point equation, 


df{x=l) I 

d\ 1-^* 


0. This yields two solutions, A±: 


A± = 


d 

2 



(4.24) 


We pick the larger one A+ as our saddle point solution, A+ = A*, as it gives a local maximum of e~^ (see Fig. |^. 
The value of e~^ diverges at A = 0, but this is not a problem as we are working in an approximation where the coils 
should have a minimum separation distance greater than Rg. If we restrict the coils to be sufficiently far apart and 
the effective radius of gyration is sufficiently small, the dominant contribution to the integral will come from A_(_. The 
solution A+ is real if d > 4i?g (where we have restored dimensional units to r). 

In the saddle point approximation, we thus find for the case of coils with distinct sequences 


{{{Fint})) RuR2 = - 


kBTalM‘^C{d,Rg)e 




V27T^6eorRj ^ 4 (^ 1 + 


where 


Cid, Rg) = 

and for the case of coils with identical sequences 




N 1 + / 


1 - 


l&Rl l&Rl 


d^ d^ 


{{{Fint))) RuR2 = - 


kBT{M^al + Mg^)C{d, Rg) e 




V2TT^{eeo)^Rl 






16 fi2 

d 2 


(4.25) 


(4.26) 


(4.27) 


The interaction energy of the two identical coils is structurally the same as that of two distinct coils, the only difference 
being that the prefactor is replaced by M'^a'l + Mg'^. A more heterogeneous polarizability sequence thus results 

in a more attractive vdW interaction between identical coils, but not between distinct coils. 


V. DISCUSSION AND CONCLUSION 

We have considered the pairwise vdW interaction for both stiff polarizable polymers and flexible Gaussian coils, and 
described their behaviors in the limiting regimes where a pair of polymers are far from each other or near each other. 
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In particular, we see that the vdW interaction between stiff polymers has a marked orientational anisotropy which 
makes it more attractive for pairs of polymers that are aligned, and moreover we find that in the near-field region, if 
the stiff polymers have identical sequences, the interaction decays as the inverse fifth power of their separation, which 
is distinct from and stronger than the inverse fourth power decay of distinct polymers. We have also seen that a more 
heterogeneous polarizability sequence also results in a more attractive vdW interaction between identical sequences, 
but not between distinct sequences, for both stiff polymers and flexible Gaussian coils. These characteristic features 
can aid the mechanism of molecular recognition between biopolymers such as dsDNA molecules in aqueous solvent. 
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